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Abstract This paper is a continuation of the previous study [Samaj, L.: J. 
Stat. Phys. 137, 1-17 (2009)], where a sequence of sum rules for the equi- 
librium charge and current density correlation functions in an infinite (bulk) 
quantum media coupled to the radiation was derived by using Rytov's fluc- 
tuational electrodynamics. Here, we extend the previous results to inhomo- 
geneous situations, in particular to the three-dimensional interface geometry 
of two joint semi-infinite media. The sum rules derived for the charge-charge 
density correlations represent a generalization of the previous ones, related to 
the interface dipole moment and to the long-ranged tail of the surface charge 
density correlation function along the interface of a conductor in contact with 
an inert (not fluctuating) dielectric wall, to two fluctuating semi-infinite me- 
dia of any kind. The charge-current and current-current sum rules obtained 
here are, to our knowledge, new. The current-current sum rules indicate a 
breaking of the directional invariance of the diagonal current-current correla- 
tions by the interface. The sum rules are expressed explicitly in the classical 
high-temperature limit (the static case) and for the jellium model (the time- 
dependent case). 
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1 Introduction 

The models studied in this paper are composed of spinless charged particles, 
classical or quantum, which are non-relativistic, i.e. they behave according 
to Schrodinger and not Dirac. On the other hand, the interaction of charged 
particles via the radiated electromagnetic (EM) field can be considered either 
non-relativistic (nonretarded) or relativistic (retarded). In the nonretarded 
regime, magnetic forces are ignored by taking the speed of light c — ► oo, so 
that the particles interact only via instantaneous Coulomb potentials. In the 
retarded regime, c is assumed finite and the particles are fully coupled to 
both electric (longitudinal) and magnetic (transverse) parts of the radiated 
field. 

One of the tasks in the equilibrium statistical mechanics of charged sys- 
tems is to determine how fluctuations of microscopic quantities like charge 
and current densities, induced electric and magnetic fields, etc., around their 
mean values are correlated in time and space. A special attention is devoted 
to the asymptotic large-distance behavior of the correlation functions and 
to the sum rules, which fix the values of certain moments of the correlation 
functions. 

Two complementary types of approaches exist in the theory of charged 
systems. The microscopic approaches, based on the explicit solution of models 
defined by their microscopic Hamiltonians, are usually restricted to the non- 
retarded regime. A series of sum rules for the charge and current correlation 
functions has been obtained for infinite (bulk), semi- infinite and fully finite 
geometries (see review [1]). The quantum sum rules are available only for the 
jcllium model of conductors (sometimes called the one-component plasma), 
i.e. the system of identically charged pointlike particles immersed in a neu- 
tralizing homogeneous background, in which there is no viscous damping 
of the long-wavelength plasma oscillations. The macroscopic approaches are 
based on the assumption of validity of macroscopic electrodynamics. Being 
essentially of mean-field type, they are expected to provide reliable results 
only for the leading terms in the asymptotic long- wavelength behavior of cor- 
relations. In general, these approaches are able to predict basic features of 
physical systems also in the retarded regime. A macroscopic theory of equi- 
librium thermal fluctuations of the EM field in quantum media, conductors 
and dielectrics, was proposed by Rytov [SH!]. 

In a recent work [5], a sequence of static or time-dependent sum rules, 
known or new, was obtained for the bulk charge and current density correla- 
tion functions in quantum media fully coupled to the radiation by using Ry- 
tov's fluctuational electrodynamics. A technique was developed to extract the 
classical and purely quantum-mechanical parts of these sum rules. The sum 
rules were critically tested on the jellium model. A comparison was made with 
microscopic approaches to systems of particles interacting through Coulomb 
forces only [Sl[7]; in contrast to microscopic results, the current-current den- 
sity correlation function was found to be integrable in space, in both classical 
and quantum cases. 

This paper is a continuation of the previous study [5]. It aims at gen- 
eralizing the previous sum rules to inhomogeneous situations, in particular 
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Fig. 1 Two semi-infinite media characterized by dielectric functions ei(u) and 
£2 (w) 



to the interface geometry of two semi-infinite media with different dielectric 
functions pictured in Fig. 1. It should be emphasized that this is not exactly 
the configuration considered in some previous studies. The standard config- 
uration was a conductor in contact with an "inert" (not fluctuating) wall of 
the static dielectric constant ew- The presence of a dielectric wall is reflected 
itself only via the introduction of charge images; the microscopic quantities 
inside the inert wall do not fluctuate, they are simply fixed to their mean 
values. Such a mathematical model can provide a deformed description of 
real materials and, as is shown in this paper, it really does. The only ex- 
ception from the described inert-wall systems is represented by the specific 
(two-dimensional) two-densities jellium, i.e. the interface model of two joint 
semi-infinite jelliums with different mean particle densities, treated in [8][9l 
ITU] . It stands to reason that in the case of the vacuum (ew = 1) plain hard 
wall, there is no charge which could fluctuate and the inert-wall model is 
therefore adequate. 

To our knowledge, the sum rules for a (fluctuating) conductor medium in 
contact with a dielectric (inert) wall obtained up to now were restricted to 
the charge-charge density correlation functions. The inhomogeneous charge- 
charge sum rules are either of dipole type or they are related to the long- 
ranged decay of the surface charge correlation function along the interface. 

The classical dipole sum rule for the static charge-charge density correla- 
tions follows directly from the Carnie and Chan generalization to nonuniform 
fluids of the second-moment Stillinger-Lovett condition [TT|I12] . The time- 
dependent classical dipole sum rule was derived in [13]. A time-dependent 
generalization of the Carnie-Chan rule to the quantum (nonretarded) jel- 
lium and the consequent derivation of the quantum dipole sum rule for the 
time-dependent charge-charge density correlations were accomplished in Ref. 

The bulk charge correlation functions exhibit a short-ranged, usually ex- 
ponential, decay in classical conductors due to the screening. On the other 
hand, for a semi-infinite conductor in contact with a vacuum or (inert) di- 
electric wall, the correlation functions of the surface charge density on the 
conductor decay as the inverse cube of the distance at asymptotically large 
distances [T5]. In the classical static case, this long-range phenomenon has 
been obtained microscopically [16, 17 as well as by using simple macroscopic 
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argument based on the electrostatic method of images [18]; the prefactor to 
the asymptotic decay was found to be universal, i.e. independent of the com- 
position of the Coulomb fluid. In the quantum case of the specific jcllium 
model, ignoring retardation effects, a nonuniversal prefactor to the asymp- 
totic decay was obtained, for both static [19] and time-dependent [141119] 
correlation functions. Recently pOpT] . by using the inhomogeneous version 
of Rytov's fluctuational theory, we have extended the quantum analysis of 
the jellium to the retarded case. We got a surprising result: for both static 
and time-dependent surface charge correlation functions, the inclusion of re- 
tardation effects causes the quantum prefactor to take its universal static 
classical form, for any temperature. The restoration of the classical prefactor 
by retardation effects was observed subsequently for arbitrary (conductor, 
dielectric, vacuum) configurations of two semi- infinite quantum media [22] . 

In this paper, we apply the inhomogeneous version of the Rytov fluctu- 
ational theory to extend the bulk sum rules, derived in [5], to the geometry 
of two joint semi-infinite media with distinct dielectric functions in Fig. 1. 
The sum rules derived for the charge-charge density correlations represent a 
generalization of the previous (dipole moment and surface charge) ones, valid 
for a conductor system in contact with an inert dielectric wall, to two fluctu- 
ating semi-infinite media of any kind. The fundamental differences between 
the results for the inert and fluctuating wall descriptions are pointed out. 
The charge-current and current-current sum rules obtained here are, to our 
knowledge, new. The current-current sum rules indicate a breaking of the di- 
rectional invariance of the diagonal current-current correlations by the inter- 
face. The sum rules are expressed explicitly in the classical high-temperature 
limit (the static case) and for the jellium model (the time-dependent case). 

The paper is organized as follows. In Sec. 2, we review the inhomogeneous 
Rytov theory of EM field fluctuations and write down basic expressions for 
the charge and current densities; explicit results for the elements of the re- 
tarded Green function tensor for the two semi-infinite media configuration 
in Fig.l are presented in Appendix. Dipole sum rules for the charge-charge 
density correlation functions are derived in Sect. 3. The sum rules related to 
the long-ranged tail of the surface charge density correlation function along 
the interface between two media are discussed in Sect. 4 which is divided into 
three parts. In the first part 4.1, we generalize the classical static analysis 
of a medium in vacuum [22j to arbitrary media configurations. Part 4.2 con- 
cerns the derivation of a classical static relation between the dipole moment 
and the large distance asymptotic of the surface charge density. Part 4.3 is a 
brief recapitulation of the quantum case, in both retarded and nonretarded 
regimes. The sum rules for the charge-current and current-current density 
correlation functions are the subject of Sects. 5 and 6, respectively. Section 
7 is the Conclusion. 

2 Fluctuational Formalism 

We consider the (3+l)-dimensional space of points, defined by Euclidean 
vectors r = (x,y,z) and time t. We shall deal with semi-infinite geome- 
tries, inhomogeneous say along the first coordinate x; it is useful to denote 
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the remaining two coordinates normal to x as R = (y,z). The model con- 
sists in two distinct semi-infinite media (conductors, dielectrics or vacuum) 
with the frequency-dependent dielectric functions e\(w) and €2(10) which are 
localized in the complementary half spaces A\ = {r = (x > 0,R)} and 
A2 = {r = (x < 0, R)}, respectively, so that the interface between the media 
is localized at x = (see Fig. 1). We shall assume that the media have no 
magnetic structure, i.e. they are not magnetoactive, and the magnetic per- 
meabilities fx% = fj>2 = 1- The two-point functions studied in this paper will 
be translationally invariant in time and in the vector space R perpendicular 
to the x axis, and so we shall use the (partial) Fourier representation 



f(t,r;t , ,T') = f(t-t',R--R , ;x,x') 



OO 



du> f dq „_i w ( t _ 4 ') +iq .( R _ R 



27T J R2 (2tt)2 



c 



)/(w,q;a!,a/),(2.1) 



where lo denotes the frequency and q = (q y , q z ) is the two-dimensional wave 
vector. 

The induced EM fields inside a material medium are random variables 
which fluctuate in time and space due to the random motion of charged 
particles. In the long-wavelength scale, the EM fluctuations are described by 
the Rytov theory [2,3,4. This theory is usually formulated in the Weyl gauge 
with the scalar potential <f>(t, r) = 0, so that the (classical) vector potential 
A(t, r) with components Aj(t,v) (j — x,y,z) determines the microscopic 
electric and magnetic fields as follows 

E= — , B = curlA (2.2) 

c at 

(we use Gaussian units), where c is the speed of light. In the context of the 
quantized EM field theory, the crucial role is played by the retarded photon 
Green function tensor D defined by 

f <^(i,r)A fe (i',r') - A k (t',r')Aj(t,r)), t > t', 
iD jk (t-t';r,r') = { (2.3) 
[0, t < f, 

where Aj(t, r) denotes the vector-potential operator in the Heisenberg picture 
and the angular brackets represent the equilibrium averaging at temperature 
T, or the inverse temperature — 1/(/cbT). For non-magnetoactive media, 
the Green function tensor possesses the symmetry 

D jk (t ~ t'; r, r') = D kj (t - t'; r', r). (2.4) 

In the Fourier space, the symmetry is expressible as 

D jk (uj,q;x,x') = D kj (uj, -q; x' , x). (2.5) 

The validity of macroscopic Maxwell's equations for the mean values of 
the EM fields implies, in the frequency Fourier space, a set of differential 
equations of dyadic type fulfilled by the Green function tensor: 

' h r ' ^ ' Afc(w;r,r0 = -47rft^tf(r-r'). (2.6) 



E 

1 



dxjdxi c 2 
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Here, in order to simplify the notation, the vector r = (x, y, z) is repre- 
sented as (x±, X2, X3). The source point r' and the index k only act as some 
fixed parameters, the boundary conditions are with respect to the field point 
r. There is an obvious boundary condition of regularity Dj k (oj;r,r') — ► 
for asymptotically large distances |r — r'| — > 00. At an interface between 
two different media, the boundary conditions correspond to the macroscopic 
requirement that the tangential components of the fields E and H = B, con- 
sidered in the gauge (|2.2[) . be continuous. The Green function tensor for the 
geometry pictured in Fig. 1 was obtained in a number of papers, see e.g. the 
method using vector wave functions [23] or the Weyl expansion method [24j 
125] . The results are usually written in a complicated way, by using the dyadic 
notation for the tensors. In order to enable the reader to reproduce easily 
calculations performed in this work, in the Appendix we present explicitly 
the Fourier transforms (|2.1|) of the tensor elements Dj k (ui, q; x, x') for two 
possible cases: the points r and r' are in the same half space or they are in 
different half spaces. 

Applying the fluctuation-dissipation theorem and assuming the symmetry 
(|2.5p , the fluctuations of the vector potential are described by the formula 

(A j (cu,( l ;x)A k (-uj,-( l ;x')) s = - coth (^yj lmD jk (u,ci;x,x'), (2.7) 

where Im means the imaginary part and (Aj(uj, q; x)A k {— u>, — q; x')) s is the 
Fourier transform of the symmetrized correlation function 

{A 3 (t,r)A k (t',r')) s = l(Aj(t, r)A k (t>, r') + A k (t>, v')A 3 {t, r)) T . (2.8) 

Here, (• • -) T represents a truncated equilibrium average, (AB) T = (AB) — 
(A)(B). 

The relation (|2.7p enables us to calculate the symmetrized two-point cor- 
relation function of arbitrary statistical quantities. Let a scalar quantity u be 
expressible in terms of the components of the vector potential, in the clas- 
sical format and in the gauge (|2.2p , as u(t, r) = J2j UjA;(t, r), where {Uj} 
(j = x,y,z) are operators acting on time and space variables. Within the 
spectral representation with a single frequency ui and two-dimensional vec- 
tor q, f(t,r) — e _lwt+lq ' R /(w, q; x), this relation takes a partially algebraic 
form u(w,q;x) = Uj \{u>, q; x)A } \{u>, q; x) (Uj can still act as operators 
on the x coordinate). It follows from the definition (|2.ip that the Fourier 
transform of the symmetrized two-point correlation function of statistical 
quantities u and v, (u(t, r)v(t', r')) s , is then determined by 

(u(x)v(x')) s w>tl = (u(w,q,a;)'u(-w,-q,x')} s 

= ^U j (u,,q;x)V k (-L,,-q;,x')(A j (x)A k (x'))Z j(l , (2.9) 
jk 

where (A j {x)A k (x')) s LJ ^ = (Aj(w, q; x)A k {-uj, -q; x')) s is given by $F§. 
The statistical quantities of our interest are the volume charge density p 
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and the electric current density j. The charge density, defined by 47rp(i, r) = 
divE, is expressible in terms of the vector potential components as follows 



p(w, q; x) = - q y A y - q z A ; _ 



(2.10) 



where the abbreviated notation Aj = Aj(u), q; x) is used. The vector compo- 
nents of the electric current density, defined by 47rj(t, r) = c curlB — <9tE, 
are expressible as 



j x (uj,q;x) 
j y (w,q; x) 
j z (w,q;x) 



c 

4tt 

c 

47T 

c 

47T 



A, 



{ i ( qyAy + qzA *> 



d a ( 2 u 

• d A 

-iq z —A x - q y q z A y 



dx 2 

2 



Ay — q y q z A 



(2.11) 
(2.12) 
(2.13) 



3 Dipole sum rules 

In the present and subsequent sections, we treat the symmetrized charge- 
charge density correlation function (p(t, r)p(0, r')) s , where we set t' — for 
simplicity. This function fulfills the obvious neutrality condition 

dr(p(t, r)p(0, r')) s = / dv'(p(t, r)p(0, r')) s = 0. (3.1) 

Let us consider a (partial) dipole moment carried by the charge-charge 
density correlation function {pit, r)p(0, r')) s , with the point r being con- 
strained to the region A\: 

/OO /"OO /> 

dx' / dx dRx(p{t,r) P (0,r')) s . (3.2) 
-oo JO J 

Note that, due to the translational invariance of the correlation function in 
the space perpendicular to the x axis, the integration J dR can be equiva- 
lently rewritten as J dR' or J d(R — R'). Interchanging naively the order of 
integrations over x' and x in (|3.2[) , rewriting J dR as J dR' and then apply- 
ing the neutrality condition (|3.1|) , the quantity Z)W(t) seems to vanish. This 
is not true. As positive x and x' become large, (p(t, r)p(0, r')) s tends to the 
bulk function (i, r — r') corresponding to the medium 1. This correlation 
function is not small when the points r and r' are close to one another. Con- 
sequently, the function in (|3.2p is not absolutely integrable which prevents 
from permuting the integrations over x' and x. Subtracting and adding the 
bulk correlation function in (13.21) leads to 



/OO />OC /• 

dx' dx dRx (p(t,r)p(0,r')) s - S 6 (1) (t, r - r') 
-oo ^0 J 

/oo />oo /> 

dx' dx / dRxSl 1] (t,r-r'). (3.3) 
-oo JO J 
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We assume that the convergence of the charge-charge density correlation 
function to the bulk function occurs on a microscopic scale, so that 



/OO pOO P 

dx' / dx / dRx\(p{t 7 r)p(0,r')) s - S 6 (1) (i, r - r')| < oo; 
-oo JO J 



(3.4) 



negative values of x' do not represent any complication for x > since the 
charge-charge density correlation function is expected to be short ranged 
along the normal to the interface. Under condition (|3.4[) . we can permute 
the x' and x integrals in the first term on the r.h.s. of (|3.3j) : regarding the 
neutrality sum rule (|3.ip . this term becomes equal to 0. Using the translation 



plus rotation invariance of the bulk correlation function S^(t, r — r') and 
the method of integration by parts, the second term on the r.h.s. of (|3.3|) can 
be easily reexpressed as follows 

/OO /"OO P -i /"OO P 

dx' / dx dRzSf'^r-r') = - / dx dRa; 2 ^ 1 ^, r) 
-oo J J ^ J -oo J 

= \f drr^W^r). (3.5) 
The calculations of this paragraph can be summarized by the equality 

/oo poo p -i p 

dx' J dx J dHx(p(t,r) P (0,r')) s = -J drr 2 ^ 1 ^r). (3.6) 

Note that this dipole sum rule depends on the bulk characteristics of the only 
one from the two media. 

We can treat similarly the (partial) dipole moment carried by the charge- 
charge density correlation function (p(t, r)/?(0, r')) s , when the point r is con- 
strained to the region A2: 



/OO p0 p 

dx' / dx dRx(p(t,r)p(0,r')y 
-OO J— OO J 



(3.7) 



The procedure analogous to the one outlined in the previous paragraph re- 
sults in 

J dx' J dx J dHx(p(t,r)p(0,r')) s = ~ J drr 2 sl 2 \t,r), (3.8) 

(2) 

where (t, r) is the bulk charge-charge density correlation function cor- 
responding to the medium 2. Combining relations (|3.6p and (|3.8p . the total 
dipole moment reads 

dx' f drx(p(t,r)p(0,r')) s = f dx' f dr (x - x')(p(t, r)p(0, r')) s 



dr r 



sW(t,r)-Sf (t,r)l. (3.9) 



We see that the dipole sum rules for an inhomogeneous configuration 
of two semi-infinite media are expressible in terms of the second moments 
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of the symmetrized charge-charge density correlation function in an infinite 
medium with the frequency-dependent dielectric function e±(u>) or e^oS)- This 
subject was studied in Sect. 3 of the previous paper [5]. The final result for the 
second-moment condition, derived by using the Rytov fluctuational theory, 
reads 



jdrr 2 ^(i,r) = j^expH 



e a (u)) ' 



(3.10) 



where the index a = 1, 2 denotes the medium. The introduced function 



(iflUJ 



coth 



(3huj 



(3.11) 



fulfills g(u>) > 1, the equality <7(w) = 1 takes place in the classical limit 
f3fioj — y 0. The integral over u> on the r.h.s. of (|3 . 1 0|) is expressible in terms of 
elementary functions perhaps only for the (one-component) jcllium model of 
conductors, i.e. the system of identical particles with the number density n, 
charge e and mass m, immersed in a neutralizing homogeneous background. 
The dielectric function of the jcllium is adequately described, in the long- 
wavelength limit q — > 0, by the Drude formula with the dissipation constant 
taken as positive infinitesimal [26j . 



e(w) = 1 - 



+ iry) 



(3.12) 



where the plasma frequency ui p is defined by to 2 = Anne 2 /m. Inserting the 
representation (|3 . 1 2|) into (|3 . 1 0[) and using the Weierstrass theorem 



lim — - — = V ( - ) =F mS(x) 
»)-»o+ x ±irj \x J 

(P denotes the Cauchy principal value), we arrive at 

r , ,„ , i 



■ J dvr 2 S b (t,r) 



2tt 



g(uj p ) cos(ujpt). 



(3.13) 



(3.14) 



In the static t = case, for all media, using path integration techniques 
and the general properties of dielectric functions in the complex frequency 
upper half-plane, the integral over ui on the r.h.s. of (|3.10[) can be formally 
expressed as [5] 



1 



- 1 



1 OO 
3=1 



1 



- 1 



Here, 



(3.15) 



(3.16) 
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are the (real) Matsubara frequencies. For the general medium composed of 
species (electrons and ions) a with the number density n a , charge e a and 
mass m a , the dielectric function fulfills the asymptotic relation ;26,27J 



e( W ) 



m„ 



(3.17) 



In the high-temperature (classical) limit (3Tilo p — > 0, each of the Matsubara 
frequencies is much larger than u p , the corresponding terms in the 

sum on the r.h.s. of (|3 . 1 5[) vanish and so the formula (|3.15[) represents the split 
of the bulk second-moment condition onto its classical and purely quantum- 
mechanical parts. We conclude that the dipole sum rules (|3.6p and (|3.8p take 
in the classical limit the following forms 



(3 



jT^d,/dB^)p(o>s--i(^- 



(3.18) 



(3.19) 



In the above derivation of dipole sum rules, only the results of the bulk 
version of the Rytov fluctuational theory were needed at the final stage of the 
analysis. In what follows we shall show how Rytov's theory can be adopted 
from the beginning in its inhomogeneous version; the true value of this ap- 
proach will be justified later. Let the point r be in the region A\, i.e. x > 0, 
the position of the point r' is arbitrary. Using the formalism of Sect. 2 and the 
explicit results for the retarded Green function tensor in the Appendix, the 
Fourier transform of the charge-charge density correlation function is found 
to be 



ei(w) 



if 



dxdx' 



5{x-x'). (3.20) 



Here, the delta function has to be understood in a macroscopic sense (dis- 
regarding microscopic structure at small distances). Taking in (|3.20[) q = 
and performing the inverse Fourier transform in time, we obtain 



/I poo 1 
dR(^,r)p(0,r')) s = - J ^ — 



i ■ \9( Ll >)t 
exp(— lut) 1m 

TTLO 

d 2 
dxdx' 



1 



(3.21) 



Since the integration by parts implies 



dx' I dx x 
'o 



dxdx 



7 5{x-x') 



= 1, 



(3.22) 



we recover the dipole sum rule (|3.6|) with the inserted bulk second-moment 
condition ()3.10j) . The dipole sum rule (|3.8|) can be verified analogously. 
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4 Surface charge density correlations 

4.1 Classical limit 

We extend the classical result of [22], Sect. II, to the configuration in Fig. 1, 
i.e., the two half-spaces Ai (x > 0) and A 2 (x < 0) filled with media charac- 
terized by static dielectric constants e\ = £i(0) and 62 = £2(0), respectively. 
We recall that e(0) — > ioo for conductors, e(0) = 1 for vacuum and e(0) > 1 
(finite) for dielectrics. We shall consider the static two-point correlation func- 
tions with zero time difference t = t', so the time variables will be omitted 
in the notation. 

For an arbitrary configuration of two points r and r' in the media, we 
shall compute the correlation function ((/>(r)</>(r')) T , where <p(r) is the micro- 
scopic electric potential created by the media at point r. It is related to the 
microscopic charge density p(r") by 

where the integral is over the whole space. In particular, we shall first cal- 
culate microscopically the electric potential due to an infinitesimal charge 
Q placed in one of the two media and then complete this calculation with 
the phenomenological electrostatics result (the method of images) for that 
potential. 

Let us introduce a test infinitesimal pointlike charge Q at point r. The 
microscopic formula for the total potential (/>tot induced at point r' is 

(<Mr')>Q = p^ + <#-')>Q, (4.2) 

where (• ■ -}q means an equilibrium average in the presence of charge Q. The 
additional Hamiltonian is H = Qtp(r). We now use the linear response theory 
for (cj)(r'))Q which says that 

(fr')) Q H^'))-W(r')ff) T 

H^r'))-/W(r>(r)) T , (4.3) 

where (■ • ■) = (• • -)q=o means the standard equilibrium average (i.e., in the 
absence of the test charge Q). Combining (|4.2[) and (|4.3[) we arrive at 

/3Q(0(r>(r)) T = [(<Mr')>Q - W)}} • (4.4) 

Now, let the point r be in region A\. According to phenomenological 
electrostatics [26], the shift of the potential average due to Q is in A\ 

where r* — (—x, R) is the position of the image charge Q' . We would like to 
emphasize that the relation (|4.5p is valid for macroscopic distances |r' — r| 
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which are much larger than the microscopic scale defined by the particle 
correlation length. If the test charge is in region A2 at r, the average potential 
{4>{ v '))q m ^2 is given by (|4.5[) with indices 1 and 2 interchanged. Finally, 
if the test charge is in region A\ at r, the average potential in region A2 is 
given by 



(&ot(r')>Q- W)) 



e 2 r' - r 



Q" 



2e 2 



ei + (2 



-Q. 



(4.6) 



A similar relation is also valid if the test charge is in A2 for the average 
potential in region A\. 

Using and (|4.5p . we obtain 



1 



1 



1 



£2 - ei 



1 



if > 0. (4.7) 



ei / |r' - r| £2 + ei ei |r' — r* 
If x, x' < 0, 1 and 2 should be interchanged. Similarly, (|4.4[) and (14.61) imply 



'\\T 



/3(0(r)0(r')) 



1 



ei + £2 / r' - r 



I I — if x > 0, x' < or x < 0, x' > 0. 

(4.8) 

The surface charge cr(R) on the plane x = at the point (0, R) is related 
to the discontinuity of the normal ^-component of the microscopic electric 
field E on the interface: 



47r<7(R)=£+(R)-£;-(R) 



(4.9) 



where the superscript + (— ) means approaching the surface through the limit 
x ^ + (x — > 0~). The surface charge correlation thus is 



'\\T 



MR)a(R')) 



JL(£+(R)£+(R') + E-(K)E-(R') 



-2 J B+(R) J B-(R')) T . (4.10) 
The electric field is related to the potential by E(r) = — V^(r), so that 



o 2 



'\\T 



(E x (r)E x (v'))' =_(0(r)0(r')> 



(4.11) 



Using (14. 7|) . we obtain 



/3(£;+(R)£;+(R')) T = 



dxdx' 



1-1 



£2 - ei 1 
e 2 + ei ei|r' - r* 



(4.12) 



Since 

d 2 1 
dxdx' |r' — 



|R-R'I 3: 



a 2 

dxdx' |r' 



IR-R'I 3 ' 
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we find 

/3(£ , ~(R)_E~(R')} T is obtained from (|4.13[) by interchanging 1 and 2. Finally, 

«(R)^(R')> T = (l - ^) j^p. (4.14) 
Using these relations in (|4. 10[) gives the classical result 

/3<a(RMR')>£ = h cl (0) = -^s (- + -- , 

|R — R | d 8tt z \ei £2 ei+£2/ 

(4.15) 

where the argument of the prefactor h c \(t — t') equals to for the considered 
static case t = t'. We recall that this classical static result is valid for asymp- 
totic distances |R — R'| much larger than any microscopic length scale. If 
in A2 there is vacuum (e2 = 1), one retrieves equation (20) in [22]. If fur- 
thermore in A\ there is a conductor (ei = 00), one retrieves the old result of 

The surface charge density a has to be understood as being the volume 
charge density p integrated along the x axis on some microscopic distance 
within the interface region. From this point of view, the formula (|4.15[) im- 
plies a sum rule for the volume charge-charge density correlation function. 
In particular, if one assumes an asymptotic behavior of the form 

0WW))l= |r!r/| 3 . IR-R'Hoo, (4.16) 
h(x,x f ) obeys, in the classical limit, the sum rule 

/oo poo 
dx' dxh(x,x') = h d (0). (4.17) 
-00 J — 00 

The above formalism can be extended straightforwardly to other geome- 
tries of the interface between media, e.g. cylindrically or spherically layered 
media, or to planarly multi-layered media. The only modification consists in 
the application of the corresponding variant of the method of images. 

It is instructive to compare the present classical result (|4.15[) . valid for 
two fluctuating media, with the previous result [IT] valid for a fluctuating 
medium in contact with the inert wall which "produces" the images, but 
does not fluctuate. For the special case of a Coulomb conductor (ei — > 00) 
in contact with the fluctuating wall of the static dielectric constant £2 = £\V: 
the prefactor h c \(0) in the formula (|4. 15[1 takes the form 

fluctuating wall: h c i{0) = — — ^ — . (4-18) 

on ew 
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On the other hand, for a Coulomb conductor in contact with the inert wall 
of the static dielectric constant ew, the prefactor /i c i(0) was found to be [17] 

inert wall: h c \(0) = --^ew- (4-19) 

We see that the two results (|4.18j) and (|4. 19|) coincide with one anther only 
for the vacuum (plain hard) wall; in vacuum, there are no charges, so that 
the description by fluctuating and inert walls should lead to the same result. 
Increasing ew beyond 1, our formula (|4.18[) predicts the suppression of the 
surface charge fluctuations while (|4. 19|) predicts their enhancement. 



4.2 Classical Surface Charge Correlations and Dipole Moment 

In the classical limit, there exists a direct relation between the dipole mo- 
ments (I3.18P , (|3.19|) and the asymptotic behavior of the surface charge density 
correlations (|4.15p . The aim of the present part is to derive this relation. 

Let us consider the potential-potential correlation function, given by (|4.7j) 
or (|4.8p . when the point r is localized at the interface, say r = 0, the position 
of the point r' is arbitrary: 

f>mm) T = (i - ^) ^- (4-20) 

Applying the Laplacian to both sides of this equation and using the Poisson 
equation A r <(j){r') = — 47rp(r'), we get 

mo)p(r')) T =(l~ ) Sir'). (4.21) 

With regard to the definition of the microscopic potential (|4.1j) . using in 
(I4.2ip the (partial, two-dimensional) Fourier transform of the Coulomb po- 
tential 

1_ I _fq_ lq . R 27r | X | 

|r| J (2^) 2e q 6 ' q ~ |q| ' [ ' 

and the convolution theorem, we get 

P |" dx {p(x)p(x'))le-^ = ± - — H— ) S(x'). (4.23) 

This equation is valid for large distances |R — R'| or, cquivalently, small 
q. Performing the small-g expansion in (14.23)) and then integrating over all 
x' G (—oo, oo), we arrive at 



/oo poo / O \ 

^ dx' J ^ dx ( P {x)p{x'))l = ±(l- J^-J 



+qf3 I dx' I dx\x\(p(x)p(x'))l =0 . 

(4.24) 
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This is the wanted relation. Inserting here the relations for the dipole mo- 
ments (|3.18j) and (|3. 19|) . we end up with 

(3 dx> j°° dx { P {x)p{x'))l = ±(1 + 1- -J—) . (4.25) 
J-oo J-oo 4?r Vei £2 ei+e 2 / 

Since, in the sense of distributions, the two-dimensional Fourier transform of 
1/R 3 is —2irq, the result (|4.25[) is equivalent to the previous one described 
by equations (f4T5|) - (l4~T7l) . 



4.3 A Short Recapitulation of the Quantum Case 

The long-range decay of the quantum surface charge density correlation func- 
tions, in both retarded and nonretarded regimes, was the subject of Refs. 
|2"TI|2"2"] . By using the Rytov formalism for a plane between two media, the 
two-point electric field correlations were derived for any point positions in 
medium 1 and 2 and the discontinuity of the electric field across the inter- 
face was related to the surface charge density. The consequent integrals over 
the frequency were treated using path integration techniques and the general 
properties of dielectric functions in the complex frequency upper half-plane. 

In the static t = t' case, the final result for the Fourier transform of the 
quantum surface charge density correlation function reads 

^U=-^(- + --—^—)+F vl (0 ) q) ) (4.26) 
47r \a e 2 ei + £2/ 

where the explicit form of the (static) function F qu (0,q) depends on the 
considered, retarded or nonretarded, regime. In the retarded regime, we have 

F Wfn ) = ^V - MiQ)-^)] 2 

qu[,q> 27r^«i(i^)e a (i6-)+«2(^i(^) eifotete) ' 

(4.27) 

where the Matsubara frequencies £j (j = 1,2,...) are defined in (|3.16p and 
the inverse lengths K a (u>, q) for the regions a = 1, 2 in (|A.1[) . For the purely 
imaginary values of the frequencies u> = i£j , the values of the dielectric func- 
tions ei,2(i£j)> an d consequently of the inverse lengths Ki^i^j), are real pos- 
itive. In the considered limit q — > 0, «i .2(1^) = £j£l,2(i£j')- Since £j cx j and, 
according to (|3.17p , £1,2 — 1 = 0(l/j 2 ), the sum in (|4.27|) converges. 

This means that the function F^} (0, q), being of the order 0(q 2 ), becomes 
negligible in comparison with the first term in (|4.26[) in the limit q — > 0. The 
prefactor associated with the asymptotic decay (|4.15|) thus reads 

^(0) = + — ^— ) • (4-28) 

87T Z \ei e 2 £1 + £2 / 

This expression, which does not depend on the temperature and the Planck 
constant, coincides with the classical result (|4. 15|) . In other words, the con- 
sideration of retardation effects makes the quantum mechanics equivalent to 
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its classical limit. The situation is fundamentally different in the nonretarded 
regime. In the limit c — > oo, it holds n a — q. We thus get from the retarded 
representation (|4.2T|) that 



i=i 



ei(i&) e 2 (iCj) ei(iCj) +ei(i^) 



(4.29) 



With regard to the asymptotic behavior ei,2(i£j) — 1 = 0(1/ j 2 ), the sum 
in (|4.29|) converges, the nonretarded function FqS^O, q) is of the order 0(q) 
and therefore contributes to the surface charge density correlation (|4.26|) . 
The prefactor hq^\o) is a complicated function of temperature. It was shown 
in [20] that for distances A ~ 1/q much smaller than c/uj p the retardation 
effects are negligible and so the nonretarded result (|4.29p takes place, while 
for A 3> c/uip the retardation results describe adequately the decay of the 
surface charge density correlations. 

In the retarded regime, the time difference between points has no effect 
on the form of the asymptotic behavior (|4.15[) . i. e. 

h$(t) = ~ (- + - - -^—) . (4.30) 



5 Charge- current density correlations 

In this section, we shall deal with the charge-current density correlation func- 
tions (p(t, r)jk(0, r')) s , where the component index k equals x = x\, y = x% 
or z = X3. We recall that for the bulk medium with the dielectric function 
e(w), these correlations were shown to satisfy the following sum rules [S] 

(3 f dr(p(t,r)j k (0,r')) s b = 0, (5.1) 

P [ drx l (p(t,v)j k (0,r , )y b =6 kl f°° ^ eX p(-i^)^Im^-. (5.2) 
J J-oc 2tt 2tti e(w) 



The static t = version of the sum rule (15. 2p is trivial for any medium: 
Since g(w)Ime _1 (ti;) is an odd function of u>, the r.h.s. of (|5.2[) vanishes. 
In the special case of the jellium model with the dielectric function (|3.12p . 
the Weierstrass theorem (13.131) permits us to express explicitly the time- 
dependent sum rule (|5.2p as follows 

J dr^(p(i,r) Jfc (0,r'))§ = 5 M ^^sm(u p t). (5.3) 

We now consider the inhomogeneous situation pictured in Fig 1. Let the 
point r' be in the region Ai (x' > 0), the position of the point r is arbi- 
trary. We start with the inhomogeneous Rytov theory (see Sect. 2 and the 
Appendix), whose results for the charge-current density correlation function 
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in the Fourier space, up to terms linear in q y and q z , can be summarized as 
follows 



(3(p(x)j x (x')) 



f\\S 



P(p(x) Jy (x'))^ 



P(p(x)j z (x')) 



9H 
27ri 

' 2tt 



Im 
Im 
Im 



1 



ei(uj) 
1 



<9.r 



6(x - x') + 0(ql,ql, q y q z ),(5A) 



ei(w) 
1 



q y 6{x - x') + 0(q y ,q z7 q y q z ), (5.5) 
q z 5(x - x') + 0(ql,ql,q v q z ). (5.6) 



Taking q = in equations (|5.4[) - (|5.6[) and regarding that 



d 



d 



dx—S(x — x') = — / dx-r— -S(x — x') = 0, 
ox 1 A ™' 



dx' 



we obtain in the leading order 

13 { dr(p(t,r)j fc (0,r')) s =0 for all k = i, y, 2. 



(5.7) 



(5.8) 



This is the analog of the bulk sum rule (|5.1[) which holds also for the point 
r' being situated inside the region A 2 . 

With respect to the equality (obtained with the aid of the integration by 
parts) 

d 

dxx—S(x — x) = —l (5-9) 
ox 

for (|5.4p and in the next order in q of e _lq ( R ~ R ) for the relations (|5.5p . 
(15.61). we find 



10) 



/? / dr^p^r^Oy)) 5 = /°° ^exp(-i^)^Im^-. (5. 

This is the analog of the bulk sum rule (|5.2|) for the point r' G A\. When 
r' S A2, an analogous sum rule is obtained by substituting in (|5. 10|1 ei (lu) by 

There exist another sum rules for the inhomogeneous situation which 
have no obvious counterpart in the bulk case. These sum rules follow from 
the application of the equalities 

/oo poo o poo poo Q 

dx dx'—S(x-x') = - dx dx'—5(x-x') = l, (5.11) 
-00 Jo ox J_ oc J dx' 

r°° r° d r°° r° d 

/ dx dx'—5(x-x') = - dx dx'—5(x-x') = -l, (5.12) 

J— oo J — 00 "X J — 00 J — oo 

to the formula (|5.4[) with q = 0. Namely, we have 

dw„_ iwt .9H im 1 



//• 00 poo pOO J 

dR / dx / dx'(p(t,r)j x (0,r')y = - — e 
JO J-00 27r 



27ri ei(^) 
(5.13) 
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and, similarly, 

/dR r dx f dx'(p(t,r)j x (0,T')r = r %-^ 9 -^-lm-^- 

(5-14) 

These relations can be verified independently by using the method for the 
dipole sum rules developed in Sect. 3. We first subtract from and add to the 
correlation function (p(t, r)jk(Q, r')) s on the l.h.s. of equations (|5.13|) and 
(|5. 14|) its bulk counterparts, corresponding to medium 1 if x' > and to 
medium 2 if x' < 0. As before, assuming that 

/>OC />00 

dR / dx dx'\(p(t,r)j x (0,r')) s - (p(t,r) Jx (0,r')) s h {1) \ < oo (5.15) 

J-oo JO 

and, similarly, 

roc />0 

dR / dx dx' \(p(t, r) Jx (0,r')) s - (p(t, r)j x (0, r'))f } I < oo, 



(5.16) 

the permutation of the x and x' integrations nullifies the contribution of the 
correlation function minus its bulk counterpart due to the sum rule (|5.8[) . Us- 
ing the translational plus rotational invariance of the bulk correlation func- 
tion in the nonzero term, we obtain 



oc 



dR / dx I dx'(p(t,r)j x (0,r'))f ) = - / dr x(p(t, r)j B (0, r'))f } 

(5.17) 



and, similarly, 



dR / dx dx'(p(t,r)j x (0y))l {2) = / drx(p(t,r) Jx (0y))f K 

J — DC J — CO J 

(5.18) 

In view of these relations, the inhomogeneous sum rules (|5.13|) and (|5. 14|) 
are in fact the consequences of the bulk sum rule (|5.2p for media 1 and 2, 
respectively. 



6 Current- current density correlations 

As concerns the current-current density correlations (jk(t, r)j/(0, r')) s (fc, I = 
x,y,z), for the bulk medium with the dielectric function e(w), they satisfy 
the sum rule [5] 

(3 [ drOfc&rJj, (0,T>))i = S kl ^ ^exp(-i^)^Im-i-. (6.1) 
J J-oo 2tt 2tt e{u>) 

In the case of the jellium model with the dielectric function (|3.12p . the Weier- 
strass theorem (|3 . 1 3|) implies 

J dr(7fc(*,r)j,(0y)>£ = fa 9 ^" cos(u> p t). (6.2) 
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A 



In the static t = case, the formula (|6.1|) can be formally expressed as 



f I uo 2 1 °° 



3=1 L 



£ 2 

J £ 2 -I- , i 2 



(6.3) 



where £j are the Matsubara frequencies defined by (|3 . 1 6[) . The first term on 
the r.h.s. of (|6.3|) represents the classical /Jfi-Wp — > limit, the second term is 
the purely quantum-mechanical contribution to the sum rule. 

For the studied configuration in Fig. 1, the inhomogcneous version of the 
Rytov method gives in the limit q — > 0, for distinct current indices, 

(3(jk{x)jl{x')) s u ^ = for k ± I, (6.4) 

for any positions of points r and r' in media 1 and 2. The relation (|6.4j) is 
equivalent to 



/ dr{i fe (i,r)i i (0 ) r')} s = for k ? I, 



(6.5) 



where the position of the point r' in media 1 or 2 is irrelevant. This is the 
generalization of the bulk sum rule (|6.f [) for k ^ I. 

Let the point r' be localized in the region Ax, i.e. x' > 0, the position of 
point r is arbitrary. For the q — > limit of the diagonal correlation function 
of the xx current components, the Rytov theory implies 



Integrating over x, this equation gives 



1 



S(x-x'). 



f Av( 3x {t,v) 3x {Q,v')Y = ~ l°° ^exp(-i^)^Im-^ 
J J-oo 2tt 2tt ex{u 



(6.6) 

r' e Ax. 

A similar expression can be derived when r' E A 2 . 

The inhomogeneous sum rules obtained up to now are quite trivial gen- 
eralizations of the corresponding bulk sum rule (|6.ip . This is no longer true 
for the diagonal correlation functions of the yy and zz current components. 
For r' g Ax, the inhomogeneous Rytov method implies 



(3 / dr(7 tf (r)j„(r')>i = (3 / dr( Jz (r)^(r')} 



Im 



2vr 



ei(w) 



(6.8) 



where the additional position-dependent function f(oj,x'), which does not 
exist in the bulk case, reads 



f(u,x') = [1-exC 



ei(w) - e 2 (w) 



ni(w)e 2 (w) + n 2 (uj)ex(oj) 
x exp I ni(w)x I 



(6.9) 
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with n a (uj) (a = 1, 2) defined by 

n%(u) = -c a (w), Rerc a (u;)>0. (6.10) 

The function f(oj,x') is equal to zero in three cases: the homogeneous case 
ei(w) = £2(w), the medium 1 is the trivial vacuum ei(u>) — 1 and far away 
from the boundary x' — ► oo. Since lim^^oo f(u>, x') — 0, the function f(u>, x') 
does not contribute to the classical limit of (|6.8p . but it does contribute to 
quantum- mechanical corrections. A similar expression can be derived when 
r' G A 2 . 

We conclude this section by noting that, according to the inhomogeneous 
Rytov theory, the interface between two media breaks up the directional 
invariance of the diagonal current-current correlations in the bulk. While the 
sum rule for the xx correlations (|6.7p has the form of the bulk one (|6.1[) . the 
yy and zz correlations l|6.8[) exhibit an additional dependence on the distance 
from the interface. We believe that this is not an artificial anomaly of the 
applied method, but the true phenomenon occurring in the current-current 
correlations functions. 



7 Conclusion 

In this paper, we applied the Rytov fluctuational electrodynamics to the in- 
homogeneous geometry in Fig. 1 to derive a sequence of sum rules for the 
charge-charge, charge-current and current-current density correlation func- 
tions. The validity of some of these sum rules was controlled independently 
by using methods developed previously in the context of the model of a 
fluctuating semi-infinite conductor in contact with an inert wall. 

In the realistic model considered here, both semi-infinite media in con- 
tact fluctuate. Comparing the classical static results (|4. 18|) and (|4. 19|) for 
the fluctuating and inert walls, respectively, we see that they coincide, as 
it should be, in the vacuum case ew = 1, but for en/ > 1 these results are 
fundamentally different. 

Some of the inhomogeneous sum rules represent a straightforward gen- 
eralization of their bulk counterparts. This is not the case of the current- 
current density correlation functions; the sum rules (|6.7p and (|6 . 8|) indicate a 
breaking of the directional invariance of the diagonal current-current density 
correlations by the interface. 
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SAS QUTE is acknowledged. 



Appendix 

In this Appendix, we present explicit forms of the retarded Green function 
tensor elements Dj^iu), q; x, x') for the two semi-infinite media geometry pic- 
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tured in Fig. 1. The half spaces Ai (x > 0) and A 2 (x < 0) are character- 
ized, besides the dielectric functions e a (ui) (a — 1,2), by the inverse lengths 
K a (ui, q) (a = 1, 2) defined as follows 



e a (u>), Re K a {oj, q) > 0. 



(A.l) 



Here, from the two possible solutions for each n a we choose the one with 
the positive real part in order to ensure the regularity of tensor elements 
Djk{to, q; x, x 1 ) at asymptotically large distances from the interface x — > ±oo. 
For simplification reasons, we shall omit in the notation the dependence of 
functions on the frequency ui and the wave number q. 

i) If the two points r, r' are localized in the same half-space, say r', r e A\ 
(i.e. x,x' > 0), we introduce a pair of functions 



u(x,x') = 



v(x, x') 



2TT?IC 2 

2-kTic 2 



e -«l|x-x'| + £1^2 - £2Kl c _ Kl(x+x >) 

£1^2 + £2«1 
e - Kl |x-x'| _ £1^2 - £2^1 c - Kl ( x+x ') 



e\K2 + e 2 Ki 

These functions satisfy the same type of the differential equation 
AttTic 2 



(A.2) 
(A.3) 



( d 2 2 



f 



uj 2 e\ 



■S(x — x'); f = u(x, x') or v(x,x'), (A. 4) 



and are related by 

du(x,x') dv(x,x r ) du(x,x') dv(x,x') 



dx dx' dx' 

The third function we shall need is defined by 

Anhc 2 k 2 — Ki 1 



w(x,x') = 



dx 



-Ki {x+x') 



(A.5) 



(A.6) 



Kl £lK 2 +£2Kl 

In terms of the introduced functions, the elements of the retarded Green 
function tensor are given by 



d 2 



D xx (x,x') = w , u(x,x') - ^j€iv(x,x'), 
dxdx' cr 



(A.7) 



d d 
D xy (x, x) = -iq y —u(x, x), D yx {x, x) = iq y —u(x, x), (A.8) 

the remaining xz and zx components are given by the replacement rule 

D xz (x,x') = D xy (x,x'){q y <-> q z },D zx (x,x') = D yx (x,x'){q y <-> q z }, 



Dyy(x, X ) — 



2 W 



u(x, x 1 ) + q 2 w{x, x'), 



(A.9) 



D zz (x, x') 



= D yy (x,x'){q y <-> q z } and 

D yz (x,x') = D zy (x,x') = q y q z [u(x,x') - w(x,x')]. 



(A.10) 
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ii) If the two points r, r' are localized in the different half-spaces, say 
r' e A\ and r e yl 2 (i.e. x' > and x < 0), we introduce the function 

a(x, x') = 47rft f e^-M* . (A . n) 

In terms of this function, the elements of the retarded Green function tensor 
are given by 

D xx (x,x') = -q 2 s(x,x') (A.12) 
D xy (x,x') = -iq y Kis(x,x'), D yx (x,x') = -iq y K 2 s(x,x'), (A. 13) 
D xz (x,x') = D xy (x,x'){q y <-» g z },D zx (x,x') = D yx (x,x'){q y <-> 

D yy (x, x) = [-ql + k 1 k 2 ]s(x, x), (A. 14) 

D zz (x,x') = D y y{x,x'){qy <-> ^ } and 

A/ Z (x, x') = £> zy (x, x') = q y q z s(x 1 x'). (A. 15) 
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